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Abstract 

Present work deals with the dynamics of vehicles which intentionally operate in 
the ground proximity. The dynamics in ground effect is influenced by the vehicle 
orientation with respect to the ground, since the aerodynamic force and moment 
coefficients, which in turn depend on height and angle of attack, also vary with 
the Euler angles. This feature, usually neglected in the applications, can be re- 
sponsible for sizable variations of the aircraft performance and stability. A further 
effect, caused by the sink rate, determines unsteadiness that modifies the aerody- 
namic coefficients. In this work, an analytical formulation is proposed for the force 
and moment calculation in the presence of the ground and taking the aircraft atti- 
tude and sink rate into account. The aerodynamic coefficients are firstly calculated 
for a representative vehicle and its characteristics in ground effect are investigated. 
Performance and stability characteristics are then discussed with reference to signif- 
icant equilibrium conditions, while the non-linear dynamics is studied through the 
numerical integration of the equations of motion. 

Nomenclature 
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AR aspect ratio 

b wing span 

c mean wing chord 

Cd drag coefficient 

Cl lift coefficient 

Ci, Cm, Cn aerodynamic moment coefficients in body axes 

Cx, Cy, Cz aerodynamic force coefficients in body axes 

F non- viscous aerodynamic force 

sum of the external forces 

F„ viscous aerodynamic force 

g gravity acceleration 

H {{Hiji)) third order apparent-mass tensor 

J inertia tensor 

h height 

h/c dimensionless height 

I vehicle length 

L transformation matrix from inertial to body frame 

m vehicle mass 

M {{Mij)) second order apparent-mass tensor 

n local normal unit vector 

^ ^ , p b q c r b . . 

p, q, r — , —, — dimensionless angular velocity components 

P linear momentum of the airstream 

Q aerodynamic moment with respect to e.g. 
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sum of the external moment with respect to the e.g. 

r = (x, z) aircraft e.g. position in Earth axes 

Yc.g. aircraft e.g. position in Body axes 

transformation matrix from a; to # 

S wing area 

Sy vehicle wetted surface 

T kinetic energy of the stream flow 

T Thrust force 

V {u, V, w) inertial velocity in body axes 

V velocity modulus 
W weight 

X, y, z Earth-fixed coordinates 

xb-i Ub, zb body axes coordinates 

Greek symbols 

a, (5 angle of attack and sideslip, respectively 

5e, 5ai 5r clcvator, ailerons and rudder angles, respectively 

5t — throttle level 

^'-max 

T] propeller efficiency 

</?,■»?, -0 Euler angles 

velocity potential 

dcf) .d(f) d(p d(p, ... . , . 

dv ^du' dv' dw^ potential by unit velocity 

7 flight path angle 
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p air density 

T = ^1 + A{h/c)'^ — 2h/c ground effect perturbative parameter 

{^xi ^y, K) = L (0, 0, 1)^ ground normal unit vector in body axes 
<^ {P-i 0.1 T^) angular velocity vector in body axes 
n engine power 
C turn rate 

Subscripts 

oo value calculated out of ground effect 

e.g. center of gravity 

a.c. aerodynamic center 

value calculated sit (p — i) — 

Introduction 

An important feature of ground effect is the influence of the attitude on the 
vehicle dynamics. This influence arises from the Euler angles, which has an effect 
on the aerodynamic force and moment coefficients. This is of paramount importance 
for all aircraft maneuvers that are performed at very low altitudes, in particular in 
phases such as takeoff and landing. 

The analysis of the aircraft stability and control during take-off and landing 
requires an accurate knowledge of the aerodynamic coefficients in ground effect^. 
Staufenbiel and Schlichting^ analyzed the longitudinal stability of an aircraft in 
ground effect. It was observed that variations of the longitudinal stability caused 
by the ground proximity is responsible for substantial changes in the landing tra- 
jectories, especially during the flare-out maneuvers. 
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More recentely the increased development of unmanned aerial vehicles (U.A.V.) 
has resulted in the consideration of these vehicles for missions at low altitudes^. 
Many U.A.Vs implement flight control programs to accomplish the required mission 
profiles and are equipped with plant control systems whose characteristics take into 
account the forces and moments that are developed in ground effect^'^. 

The vehicles which are considered in this study are the "wing-in-ground-effect" 
(W.I.G.) crafts, also called Ekranoplans, that are high-speed low-altitude flying ve- 
hicles that purposely utilize the favorable ground effect. W.I.G. craft are of interest 
due to the peculiarities that they exhibit in comparison to conventional airplanes. 
Within the framework of high-speed low-altitude transportation, the high aerody- 
namic efficiency that these vehicles achieve together with low wing aspect ratios and 
low structure weights make them more suitable than conventional aircraft^. 

During horizontal flight with (p — (3 — 0, the ground plane imposes a boundary 
condition on the aerodynamic field that symmetrically reduces the downward flow 
about the aircraft. This results in a downwash reduction at the tail together with 
an increase in the lift slope of both wing and tail^ that depends on the height above 
the ground. However if the vehicle flies with an arbitrary orientation with respect 
to the inertial frame, the ground alters the vehicle aerodynamics in such a way 
that the pressure distribution on the aircraft depends on both height and attitude. 
Therefore the aerodynamic coefficients, which are a function of both height and 
aerodynamic angles, also vary with the Euler angles. This effect, which modifies 
performance and flying qualities, is observed by pilots when an aircraft is flown 
close to the ground. For example, during flight tests on W.I.G. vehicles^, horizontal 
banked turns were carried out to demonstrate that the sideslip, rudder and ailerons 
angles are quite different from those needed out of ground effect. While these angles 
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are quite small for steady banked turns out of ground effect, for a W.I.G. craft 
they are quite large. Kornev and Matveev^, considered the effect of the roll angle 
and proposed semi-empirical expressions for the aerodynamic coefficients which are 
based on results from numerical simulations. They noted that, during a banked turn 
near the ground, the vehicle develops additional aerodynamic forces and moments 
that depend upon the roll angle. 

Prom a safety perspective, a banked turn must be performed with a limited roll 
angle and at a high enough altitude to avoid the risk of touchdown. Thus, a special 
maneuver strategy is adopted to maintain the distance between the wing tip and 
ground surface.^ This maneuver requires a very large turning radius that strongly 
reduces the maneuverability of the aircraft. For this reason it is prudent to use 
a control system to minimize the turning radius^. The design of such controllers 
requires a knowledge of the aerodynamic coefficients that are developed in ground 
effect^, while also accounting for the influence of the attitude on aircraft's forces and 
moments. 

Another source of the significant variations of the aerodynamic coefficients of a 
W.I.G. craft, is nonzero sink rates. The vertical velocity of the aircraft determines 
the unsteady or dynamic ground effect which, in turn produces aerodynamic forces 
and moments that depend upon the flight path angle. Such an effect occurs during 
W.I.G. maneuvers as in the case of the vertical jumping^, which is very short climbing 
phase that is used to avoid coUisions with low obstacles. 

In the present work the flight dynamics of an aircraft in ground effect are ex- 
amined through an analytical procedure that estimates the aerodynamic forces and 
moments and takes into account the aforementioned effects, such as attitude and 
sink rate. 
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In previous work Chang and Muirhead^ experimentally examined the forces on 
low aspect ratio wings with sink rate. Nuhait and Mook^ numerically examined the 
phenomena using an unsteady vortex-lattice method. Nuhait and Zedan^ analyzed 
the unsteadiness induced by the vertical velocity using a vortex lattice method, with 
a predetermined wake shape. More recently Han and Yoon^° investigated the 2D 
ground effect of flat plates in tandem configuration using a discrete vortex method. 
They observed that the unsteadiness has a marked impact on the performance of a 
tandem configuration. 

In Ref. 11, a mathematical model of flow about lifting bodies, in steady and 
unsteady ground effect is presented. The linearized equations of motion for the 
longitudinal dynamics, in which the aerodynamic coefficients are expressed by means 
of linear derivatives with respect height, pitch angle and their time derivatives are 
developed. Kumar^^, examined the dynamics of a W.I.G. craft including the effect 
of a perturbation in the forward velocity; the stability of the W.I.G. craft was 
also examined from an eigenvalue analysis of the linearized equations of motion. 
Staufenbiel^^ analyzed the longitudinal stability of aircraft in ground effect using 
the characteristic equation. The nonlinear effects were also examined by numerical 
calculation of the time-histories of the motion equations. 

Although other references to the aerodynamic forces and moments in the ground 
effect ^^^^^ can be found, to the author's knowledge the effects exerted by the atti- 
tude of the aircraft on the flight dynamics have not been examined. Therefore the 
objective of the present work is to develop an accurate mathematical model that 
may be used for the analysis of the performance and flying qualities of aircraft in 
ground effect. 

In the present work the aerodynamic forces and moments are derived from the 
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Lagrange equations. The Lagrangian function of the physical problem is given by 
the kinetic energy of the airstream, which is expressed in terms of the aerodynamic 
and Euler angles. This approach allows the forces and moments in ground effect to 
be accurately formulated and the variations in the performance and flying quality 
to be interpreted. 

Equations of Motion 

In this section the rigid aircraft equations of motion, which are used to study the 
vehicle dynamics in ground effect, are presented. These equations are written in the 
form^ 

m(v + u? X v) = 

Ju> + u> X Jw = (1) 
r = L^v 

4 = R-^w (2) 

The elements of the diagonal matrix J are 

J^^ = 0.02 m 6^ Jyy = 0.02 m J^, = 0.02 m (6^ + Z^) 

=: T+m Lg+F+F^ and = Q are the sum force and the moment with respect 
to center-of-gravity of the external forces, where m Lg represents the weight force 
expressed in body axes. F^,, which has a purely dissipative action, is the aerody- 
namic viscous force that is given as F^ = — pV^SC d / (cos a cos /3, sin /3, sin a cos /?) , 
with Cof — 0.047 that does not depend on ground effect. Thus F + F^ and Q are, 
aerodynamic force and moment which, in body axes are expressed as 
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The aircraft is assumed to be equipped by a propulsive system which is a constant 
speed propeller driven by a reciprocating engine. Thus, the thrust force T, which 
is aligned with xb, is calculated by dividing the propulsive power Umax^rV by the 
velocity component u, where the propeller efficiency rj is assumed to be constant 
and equal to 0.82. 

Analysis 

The goal of the present section is to examine how the attitude acts upon the 
flow-field about a vehicle fiying in ground effect. Consider Fig. 1, a rigid aircraft 
which moves in a potential flow and in close proximity to a rigid surface of inflnite 
extent. The attitude is represented by the vector 

A=(— sint?, sin cost?, cos cost?) (4) 

which is the normal unit vector of the ground surface in body frame. The wakes 
behind the body are considered to be rigid and their shedding lines are assigned. 
Therefore the velocity potential is the scalar product between the derivative of the 
potential with respect to the velocity and the flight velocity v in body frame^^, i.e. 

(j) — ^ ■ V = • (cos a cos /3, sin /5, sin a cos /3) V (5) 
av ov 

Since the velocity potential satisfies the Laplace equation and satisfies the integral 
equations derived from the application of the Green's second identity, the derivative 
of the potential with respect to the velocity is given by the sum of surface integrals^'''. 
This derivative can be written in terms of A as 

where ^ = ((— — ttt-)) is the matrix of the second order derivatives of the 
awaX ovioXj 

potential. 
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While Eq.(5) yields the relationship between and the aerodynamic angles, the 
second term of Eq. (6) accounts for the effect of the attitude on the vehicle aerody- 
namics. The kinetic energy of the airstream using Eqs. (5) and (6), is 

T^IpJ v-n0d5 = ivMv+ iv{H(A-Ao)} V (7) 



where 



M 



/ Is, 



ndS, H 



'sv dvdX 



ndS 



(8) 



Isv dv 

are second and the third order tensors. Eq. (7), describes the kinetic energy of the 
stream in terms of aerodynamic angles; the effect of the vehicle attitude is accounted 
for through the second term, where the quantity H(A — Aq) is a second order tensor 
that is calculated as H(A — Aq) = Hijk{Xk — A^q)- 

It is worthwithe to highlight some characteristics of the kinetic energy. First, 
for flat ground, the heading angle does not modify the aerodynamic field about the 
aircraft, therefore in Eq. (7) the kinetic energy does not depend on tjj. Furthermore, 
because of the vehicle symmetry about the plane {xb, zb), T is an even function of 

and j3. Therefore both M and H are expressed in the form 



M 








Mi3 





M22 





Mi3 





M33 



; H 



hii hi3 

h22 
hi3 h33 



(9) 



with hjj- = {hiji, 0, hijz), i,j ^ 2. 

Eqs. (8) establish that M and H are functions olh/c that do not vary with the 
Euler and aerodynamic angles. In the present work, both the tensors M and H are 
expressed in terms of the dimensionless height through the parameter r 



/ = /oo + /lT + /2T^ 

where r is defined in the nomenclature. 



(10) 
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Calculation Method 

In the equations of motion, the representation of the aerodynamic force and 
moment coefficients in terms of the Euler angles, flight path angle and dimensionless 
height, is necessary to interpret the performance and flying qualities in ground effect. 
This section deals with the calculation of aerodynamic actions developed by a vehicle 
in ground proximity. The procedure presented here allows the calculation of the 
aerodynamic forces and moments through the Lagrange equation method. This 
approach has been already applied to the calculation of the aerodynamic forces and 
moments of ultralight aircraft that are flown in the presence of a wind gradient^^. 
To derive the expression of force and moment, the kinetic energy of the potential 
flow that is generated by the vehicle is considered. Then the aerodynamic force F 
and moment Q are calculated using the Lagrange equation method in the general 
form^^'^^ 

^ ddT dT 
dT 

Q = Qa.c. + (ra.c. - rcff.) X F = -V X — + (r„.e. - V^,g) X F (12) 

Substituting into Eqs. (11) and (12) the expression of T that is given by Eq. (7), 
the aerodynamic force and moment are 

F = -[ M + H(A - Ao)]v - HAv - w x [ M + H(A - Aq) ]v (13) 

Q = -V X Mv - V X [H(A - Ao)]v + (r„.e. - r,.^.) x F (14) 

Eq. (13) states that F is the sum of three terms, the flrst of which contains the 

two quantities M and H that are 

dm dm . dMi. dMii . 
^=d^"'^~dh ^''""^ ^d^'^'^'dh' 

9H 9H . dHijk dHijk . 
H = -— v + — Vsm-f = vi + — ^ ysm7 

dr^ dh driy, dh 



(15) 
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According to the classical lifting bodies theory, M and H are proportional to the 

velocity circulation around the aircraft. They represent the time derivatives of M 

and H which are related to the shedding wake surfaces past the aircraft and to the 

smJi rate. — — and — are the apparent mass terms relative to the unit oi wake 

length which gives the aerodynamic force developed in horizontal flight, whereas 

both and produce the unsteady ground effect when 7 7^ 0. This last result 

is a consequence of the observation that the pressure on the vehicle surface depends 

V ■ V d(j) 

on the sink rate through the Bernoulli theorem p = p{const 

V is the local velocity of the stream, while ^ = ^ V sin 7 is the unsteady term 

ot oh 

caused by the sink rate. The second term in Eq. (13) depends on A, and represents 
a force that varies with the angular velocity. In fact it can be written taking into 
account that A is related to the angular velocity through Eq.(2) 



d\ ■ dX ^ , 
9$ 9# 



(16) 



where 



dX 
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dXx 
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dX, 


dX, 


dX, 



-C0S1? 

cos (p COS —simp sin ■& 
— sin ip cos ■»? — cos (p sin 1? 



(17) 



dcp d'd dijj 

Hence, this is a contribution to the rotational derivatives that tend to zero as h/c ^ 
00. Also the last term of Eq. (13) gives a contribution to the rotational derivatives, 
which in turn depends upon the height and vehicle attitude. 

In Eq. (14) for the moment, the first term is a function of a and /3, whereas the 
second term is a function of the attitude variations. 
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Thus, force and moment are 




V + [^7-l^sin7 + 



(9H 



v](A-Ao)}v 



(18) 



{H 



R-^aj}v - uj X { M + H(A - Ao)}v 



Q = -V X {M + H(A - Ao)} V + (r„.e. - r^.^.) x F 



(19) 



Therefore, Eqs. (18) and (19) allow the determination of F and Q in ground prox- 
imity and account for the attitude and sink rate. 



To validate the method, several comparisons with existing data in the literature 
are presented. 

The aerodynamic force and moment coefficients are evaluated from Eqs. (13) 
and (14), where M and H are numerically calculated using Eq. (8). The velocity 
potential is computed by an unsteady vortex-lattice code, developed by the author, 
that takes into account the ground effect using the method of images. Specifically the 
presence of the ground is simulated by placing a specular image of the lifting body 
at an equal distance below the ground plane. Thus, two symmetrically positioned 
lifting bodies are considered to determine a flat ground surface. 

The results of the first case are presented in Fig. 2a. The figure shows increments 
of the lift coefficient in terms of /i/c for a delta wing having an aspect ratio equal to 
1.456; h is measured at the midpoint of the root chord. The simulation is computed 
for an angle of attack of 22.1° and a fiight path angle equal to —10°. The dashed 
line represents the data of Ref. 8, while the symbols are the experimental data of 
Chang and Muirhead^. The figure shows that the hft coefficients obtained by the 
present method, shown by the solid line, are in good agreement with the those of 
both the Refs. 7 and 8. 



Validation of the Method 



13 



Figs. 2b, 2c, and 2d show, respectively, the variations of hft, drag and pitching 
moment coefficients for a delta wing having an aspect ratio of 1.456, in terms of the 
dimensionless height, measured at the trailing edge. The simulations are made for a 
= 10° and 7 = 0, —20°. The difference between the values calculated by the present 
method and those reported in Ref. 8 is always less than 6 %. 

The next case, shown in Figs. 2e, 2f and 2g shows the effect of aspect ratio on 
the aerodynamic coefficients. The plots show three different calculations for AR — 
1.456, 1.072 and 0.705. According to Ref. 8 the height is measured at the trailing 
edge and the solid symbols represent the data given in Ref. 8. The calculations 
show the unsteady ground effect, that is obtained for a = 10° and 7 = -20°. The 
results are in good agreement with those of Ref. 8. 

These comparisons demonstrate that the present method provides an accurate 
estimation of the aerodynamic coefficients for the examined wings in steady and 
unsteady ground effect over a wide range oi h/c variations. 

Reference Vehicle 

Consider a model of a vehicle sketched in Fig. 3, the main parameters of which are 
reported in Table 1. The vehicle is chosen according to the "Lippisch Design" that 
combines a high positioned tail with an inverted delta wing having negative dihedral 
along the leading edge. Such configuration is longitudinally stable at different h/c, 
as a result of an adequate aerodynamic moment developed by the tail which is 
located out of ground effect. 

The aerodynamic coefficients of the aircraft are estimated using Eqs. (18) and 
(19), where the matrices M and H are numerically evaluated by means of Eqs. (8). 

The velocity potential is calculated by the boundary element code VSAERO^° 
that is capable of solving the complex aerodynamic field around an aircraft in the 
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presence of the ground. 

The incremental aerodynamic coefficients caused by the control angles are as- 
sumed to be constant in any situation. They are 

AC,, = -0.029 5e^ - 0.0493 - 0.113 

/^Cy = 0.157 5r 

AC^ = -0.32 Se - 0.0112 Sr'^ + 0.1352 

(20) 

ACi = -0.07 5a + 0.04 5r 

ACm = -0.923 Se + 0.20 - 0.0055 

ACn = +0.0035 5a - 0.072 5r 
where the control angles are expressed in radian. 

The results that deal with the vehicle characteristics in ground effect are next 
presented. In the subsequent figures the lines represent the calculation made using 
Eq. (18) and (19), whereas the symbols indicate the aerodynamic coefficients derived 
from surface integrals of the pressure forces computed in the code. 

Fig. 4 shows how the aerodynamic characteristics vary with the height. The lift 
coefficient and aerodynamic efficiency in steady ground effect are shown in terms of 
the angle of attack for h/c = 0.167, 0.614, 2.792, oo, with ip — (5 — 0. The variations 
in Cl are similar to the corresponding data presented in Ref. 1, while Cl/Cd, whose 
maximum value varies from about 7 (out of ground effect) to more than 12 (h/c = 
0.167), is in good qualitative agreement with the data presented in Ref. 11. 

The influence of the pitch angle on the aerodynamic coefficients in steady ground 
effect is shown in Fig. 5 for Cl, Cd and vs. i?. The calculations are carried out 
at h/c =1, with j = j3 — ip — 0, for the two cases H = (dashed lines) and H ^ 
(solid hues). The case with H = corresponds to the usual assumption of the 
aerodynamic coefficients that do not depend on ip and -d, whereas the calculation 
with H 7^ accounts for the attitude effects. It is worthwhite to emphasize that 
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the pitch angle produces marked variations on all the aerodynamic coefficients. In 
the case H 7^ 0, Cl and Cd exhibit significant increments when > 0, while the 
attitude acts upon the pitching moment so as to result in a quasi-linear variation of 
Cmi'd)- For each coefficient, the maximum percentage differences between the two 
cases can reach 20% of the corresponding value calculated at = = 0. 

The effects of the lateral attitude are next examined. Fig. 6 shows all the 
aerodynamic coefficients in body axis in terms of roll angle, calculated for a — h° 
and /3 = 7 = 0. Because of the vehicle symmetry, the longitudinal coefficients, 
C-r, Cz and Cm and the lateral ones, Cy, Ci and C„ are, respectively, even and odd 
functions of ^p. For \ip\ < 10°, the increment in the longitudinal coefficients can be 
greater than 15% for h/c < 1. Note that the variations of Cx and correspond 
to increases oi Cl/Cd that can lead to an improvement in the vehicle performance 
during turning maneuvers. As far as the lateral coefficients are concerned, linear 
variations are observed, for h/c < 1 that can be greater than those caused by a 
sideslip angle = \(p\. Therefore, the attitude effects modify the force and moment 
coefficients; the degree of change depends upon the dimensionless height. 

The influence of the flight path angle on the aerodynamic coefficients is next 
examined. Fig. 7 shows C^, Cd and Cm in terms of h/c, for 7= -5°, 0, +5°. In 
these results, for which a = 5°, /3 = </? = 0, it is evident that the ground effect is 
stronger for negative flight path angles at relatively low h/c. 

The proposed model also includes the effects of angular velocity. Force and mo- 
ment rotational derivatives are calculated from the derivative of Eqs. (18) and (19) 
with respect to p, q and r. Fig. 8 shows the dimensionless longitudinal rotational 
derivative in terms of the angle of attack, calculated at /3 = 99 = 0, for various h/c. 
Large variations in the derivative can be seen. 
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Comparisons between the present results and those reported in Kornev^ and 
Staufenbiel^^ are next presented. In Ref. 5, the rotational derivative Ciq of the 
W.I.G. craft ELAOl is given. Figure 9 shows the comparison between the deriva- 
tives obtained in Ref. 5 (horizontal axis) and these calculated using the proposed 
method (vertical axis). In Fig. 9a the derivative Ciq calculated at h/c — 0.1 is 
represented, whereas Fig. 9b reports Cl'^ and Cm-y- Although the vehicle is geo- 
metrically different with respect to the reference vehicle, the present results are in 
relatively good agreement with the data of Ref. 5. As for Ref. 13, the data concerns 
the W.I.G. craft X-113, which is also designed according to the Lippish criteria. Al- 
though the two vehicles have some geometrical differences. Table 2 shows that the 
aerodynamic coefficients calculated by the present method are in good agreement 
with those from Ref. 13. 

Results and Discussions 

In what follows, to study the influence of the ground effect on state and control 
variables at trim, some significant situations which correspond to the level and turn- 
ing flight will be considered. Furthermore the vehicle stability will be investigated 
by means of the eigenvalues analysis applied to the linearized motion equations, 
whereas the nonlinear analysis of the vehicle motion will be made integrating the 
full set of the equations of motion. 

Eqs. (1) and (2) are used as equations of motion and the contribution of the 
three control angles to the aerodynamic forces and moments is taken into account 
through incremental force and moment coefficients as reported in Eqs. (20). 

As a first result. Fig. 10 summarizes the trim condition in level flight for /3 = 
at different h/c. The diagrams show the angle of attack, elevator angle and required 
throttle in terms of the flight speed. 
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Since in ground effect tfie lift coefficient increases wfien h/c diminisfies, tfie angle 
of attack decreases with h/c (Fig. 10a). 5e (Fig. 10b) monotonically increases with 
V , which indicates that there is a stable behavior except the lower velocities and 
altitudes, where the curves show decreasing branches {h/c < 1). This last feature 
only occurs at high angle of attack and is the result of the simultaneous effects of ■»? 
and a on the aerodynamic forces and pitching moment. Fig. 10c shows the required 
throttle as a function of the flight speed. Each curve exhibits a distinct minimum 
value that increases with h/c. Since the ground effect is stronger at high angle of 
attack, the throttle for /i/c= 0.5 and 1 is less than the throttle for /i/c — > oo only 
for V < 55 m s~^. 

In order to evaluate the attitude effect on symmetric flight situations, two different 
trim conditions in straight flight, aXh/c= 1, are examined (Fig. 11). The continuous 
and dashed hues indicate the solutions calculated for H 7^ and H = 0, respectively. 
The two solutions show substantial divergences which are caused by the different 
aerodynamic coefficients in the two cases. At lower flight speeds the angle of attack 
shows important differences between the two solutions (Fig. 11a), whereas minor 
divergences are evident in Fig. lib, where the elevator angle is displayed. The httle 
angles of attack obtained for H 7^ are in agreement with the data given in Fig. 
5 that shows relatively high lift coefficients when the attitude effects are accounted 
for. As a result, for H 7^ 0, lower flight speeds than those calculated for H = 0, are 
obtained. Also the required throttle (Fig. 11c) depends upon the attitude. In fact, 
d acts on the aerodynamic forces in such a way that the throttle exhibits significant 
reductions for velocities less than 40 m s~^. 

The third case, shown in Fig. 12, analyzes the effect of the roll angle on a banked 
maneuver. The aircraft carries out a "truly banked" turn in the horizontal plane. 
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and the angular velocity (p, q, r) = (— sini?, cos ■»? sin (/?, cos ■»? cos </?) C is vertical 
with C = 5°s~^. In this maneuver the sum of weight and centrifugal force at the 
aircraft e.g., is in the vehicle symmetry plane {xb,zb)- In order to avoid possible 
touchdown during the maneuver, it is assumed that /i/c = 2 in all situations. In 
the diagrams the solutions which correspond to H = and H are shown. Fig. 
12a shows that there are higher differences in the angles of attack between the two 
cases at lower flight velocities. As in the case of a "truly banked" maneuver that 
is performed out of ground effect, the solution for H = shows very small sideslip 
angles (Fig. 12b), whereas for H 7^ nonzero sideshp are observed. The reason of 
such discrepancy is due to (p that causes additional force and moment terms that 
are thus balanced by nonzero sideshp angles. The roll angle (Fig. 12c) exhibits 
little differences between the two solutions over the entire range of the flight speed, 
whereas more significant differences are seen in the pitch angle (Fig. 12d). If H = 0, 
ailerons and rudder are very little as in a banked turn at high altitude, whereas, for 
H 7^ the attitude effects lead to a disagreement between the two solutions. Fig. 
12h shows the effect on throttle for both the cases. As a result of the aerodynamic 
efficiency increments caused by the nonzero roll angles (see Fig. 6), lower throttle 
levels are needed for H 7^ 0. This means that, at least in these flight conditions, 
the drag increment due to the maneuver is in part counterbalanced by the drag 
reduction caused by ip. 

The vehicle stability is next evaluated through the eigenvalue analysis applied to 
the linearized equations of motion, where the various modes are identifled by means 
of the eigenvector analysis. The root locus is calculated for straight flight at the 
speed of 40 m/s, varying h/c from 00 to about 1. 

The attitude effects on the vehicle stabihty are summarized in the two sets of Figs. 
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13a, 13b and 13c, 13d which show, respectively, the two cases H = and H 0. 
Because of the different aerodynamic force and moment coefficients developed in the 
two situations, the eigenvalues vary with /i/c in a different manner. However, out of 
ground effect, the vehicle shows three oscillating stable modes (which are phugoid, 
short period and Dutch roll) and two aperiodic modes (that are the roll convergence 
(stable) and the spiral mode (unstable)). 

The case H = is first discussed. Phugoid is a stable mode in the entire range 
of /i/c, whose eigenvalues significantly vary for h/c < 4 (see Fig. 13a and 13b). 
Its imaginary and real parts, respectively, increase and diminish as /i/c — > 0. At 
each h/c the short period is a stable mode that has significant variations in the 
eigenvalues for h/c < A that are due to the nonhnearities of Cm{'&) (see Fig. 5). 
The Dutch roll eigenvalues show smaller variations. Their imaginary part do not 
change with h/c, while the real part displays more marked variations. As for the non 
oscillating modes, the positive eigenvalue associated with the spiral mode increases 
as /i/c — > 0. Although the rolling convergence remains a stable mode, it presents a 
negative eigenvalue that rises as soon as h/c diminishes. 

Figs. 13c and 13d show the root locus for H 7^ 0. The attitude effect produces 
significant differences with respect to the previous case. The phugoid mode exhibits 
the same behavior with the exception of lower h/c, where it becomes unstable (See 
Fig. 13c). Such instability is produced by the combined effect of both pitch angle 
and height upon the lift coefficient which is more pronounced at lower h/c. With 
respect to the case with H =0, the attitude infiuences the pitching moment such that 
the derivative C^^^ is almost constant with the height; therefore the short period 
eigenvalues have smaller variations with h/c. 

Pronounced discrepancies between the two cases are evident in terms of the lateral 
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modes. It is seen that the Dutch roll eigenvalues are quite different from those 
of the preceding case. They show increments for h/c < 4 that make the mode 
unstable at h/c^ 1.5. More important divergence are apparent for spiral mode and 
rolling convergence. The Spiral mode is completely changed. It becomes stable at 
about h/c — 20, whereas, for h/c ~ 10.5, it degenerates, together with the rolling 
convergence, in a new mode called oscillating rolling that is the consequence of the 
roll angle influence on the aerodynamics force and moment. The oscillating rolling, 
which occurs under a certain h/c, is a stable oscillating mode that modifies the 
lateral stabihty. Also the roUing convergence depicts nonneghgible variations. Its 
eigenvalue strongly varies for 11 < h/c < 20, while, as seen, for h/c ~ 10.5 it 
disappears. 

In the reported calculations the attitude considerably changes the lateral stability 
of the vehicle. The stabihzation of the spiral mode, the disappearance of both spiral 
and rolling convergence and the subsequent genesis of the oscillating rolling yield a 
rather different lateral stability than one where the attitude is not accounted for. 

In order to estimate the dynamic response, time histories, which are calculated 
by integrating the full set of the motion equations, for H are next examined. 
The first simulation takes place in the vertical plane at various initial heights ho, 
where the vehicle, initially in level flight at equihbrium at the speed of 40 m s~^, is 
perturbed by a pitch rate of 0.8 s~^. Fig. 14 gives the time histories of the perturbed 
height and angle of attack. In the diagrams phugoid and short period modes can 
be seen. Note that the characteristic times of both the modes are in agreement 
with the previously calculated eigenvalues. In all the cases the phugoid mode is 
persistent during the entire simulation period, whereas the short period, due to its 
strong damping, only acts at the beginning of the motion. As seen in Fig. 13, as 
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soon as Hq/c diminishes, the phugoid period decreases and its amphtude is more 
damped. 

The second case, shown in Fig. 15, concerns a perturbed lateral motion, where 
the initial condition is the same one of the previous case, while the perturbed action 
is now given by a roll rate of 0.2 s"-*^. Out of ground effect ( Hq/c — > oo), it is 
possible to identify rolling convergence, Dutch roll and spiral mode. For each initial 
height, the time- histories show characteristic times that are in agreement with the 
data obtained through the eigenvalues analysis. For /iq/c = 4, the oscillating rolling 
mode is apparent together with the Dutch roll, while both rolling convergence and 
spiral modes disappear. It is interesting to note that, at /iq/c = 2, the Dutch roll is 
obscured by the oscillating rolling which produces the main effect. This result can 
be also obtained by analyzing the eigenvectors of each mode. 

Conclusions 

In this study the dynamics in ground effect is investigated accounting for the 
influence of the attitude on the vehicle dynamics. The Lagrange equations method 
applied to lifting bodies is used to express the forces and moments that are devel- 
oped in ground effect. This method has general validity and contains the important 
elements for describing ground effect. In particular the method leads to the deter- 
mination of aerodynamic force and moment coefficients in terms of roll, pitch and 
flight path angles. Comparisons between existing data in the literature and those 
obtained with the present method show that the procedure provides accurate results 
for aerodynamic force and moment of wings in steady and unsteady ground effect. 
Once the method is validated, the aerodynamic coefficients of a vehicle geometry are 
examined. The attitude effects on the vehicle aerodynamics are identified and the 
unsteadiness induced by nonzero sink rates are explained. The trim calculation for 
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the banked turn clarifies the behavior observed in ground effect, while the stabihty 
analysis yields explanations about the flying qualities with particular reference to 
the spiral mode stabihzation and to the existence of the oscillating roUing. The non- 
linear analysis, that is obtained by integrating the full set of the equations of motion, 
gives results in accordance with those obtained by means of the linear stability. 
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Table 1. Dimensions and mass 



Overall Length, I (m) 6 

Wing Span, b (m) 6.8 

Planform Area, S (m^) 16 

Maximum Power, Umax (kw) 200 

Overall Weight, W {N) 9806 
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U 1 

h/c 


OO 




1.0 




0.4 






(3.0) 


3.47 


(3.2) 


3.7 


(3.6) 


4.02 


Clh 


(0.0) 


0.0 


(-0.035) 


-0.042 


(-0.35) 


-0.38 


Cm a 


(-0.68) 


-0.7 


(-0.70) 


-0.73 


(-0.73) 


-0.76 


Cmh 


(0.0) 


0.0 


(0.003) 


0.0037 


(0.055) 


0.063 




(-4.44) 


-4.20 


(-4.44) 


-4.03 


(-4.44) 


-3.97 



Table 2. Dimensionless aerodynamic derivatives vs. h/c. The values in the 

parenthesis are from Ref. 13 
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List of captions 



Fig. 1 Influence of </? and ■»? on a vehicle in ground effect. 

Fig. 2 Comparison of the results, 
a) Delta wing in dynamic ground effect. AR — 1.456, a — 22.1°, 7 = —10. Symbols 
are from Ref. 7. b), c), d) Delta wing in ground effect at different 7. AR — 1.456, 
a — 10°. e), f), g) Delta wings with different aspect ratios at a — 10°, 7 = —20°. 
Symbols are from Ref. 8. 

Fig. 3 The reference vehicle. 

Fig. 4 Vehicle lift coefficient and aerodynamic efficiency in ground effect obtained 
for h/c = 0.167, 0.614, 2.792, 00, at = /3 = 0. 

Fig. 5 Influence of pitch angle on the aerodynamic coefficients, h/c — 1, = P = ~ 0. 

Fig. 6 Influence of roll angle on the aerodynamic coefficients, a = 5°, /3 = 7 = 0. 

Fig. 7 Influence of flight path angle on the aerodynamic coefficients, a = 5°, /? = = 0. 
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Fig. 8 Vehicle rotational derivatives in ground effect in function of the angle of 
attack. (3 = ip = 0. 

Fig. 9 Comparison of the dimensionless aerodynamic derivatives for /i/c = 0.1. Data 
from Ref. 5 in horizontal axis, present results in vertical axis. 

Fig. 10 Trim calculation in horizontal flight at different heights. P = ip — 0. 

Fig. 11 Influence of the attitude on the trim in horizontal flight. 

Fig. 12 Trim analysis of a "truly banked" turn with C = 5° s~^ 

Fig. 13 Root Locus. H = a) and b), H c) and d). 

Fig. 14 Time history of the perturbed longitudinal motion beginning from different 
initial heights, (H 7^ 0). 

Fig. 15 Time history of the perturbed lateral motion beginning from different initial 
heights, (H^ 0). 
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Figure 1: Influence of 99 and on a veliicle in ground effect. 
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Figure 2: Comparison of the results, a) Delta wing in dynamic ground effect. 
AR = 1.456, a = 22.1°, 7 = -10. Symbols are from Ref. 7. b), c), d) Delta 
wing in ground effect at different 7. AR — 1.456, a = 10°. e), f), g) Delta wings 
with different aspect ratios at a = 10°, 7 = —20°. Symbols are from Ref. 8. 
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Figure 3: The reference vehicle. 
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Figure 4: Vehicle lift coefficient and aerodynamic efficiency in ground effect obtained 
for h/c = 0.167, 0.614, 2.792, oo, at = /3 = 0. 




Figure 5: Influence of pitcli angle on the aerodynamic coefficients. /i/c=l,7 = /3 
= = 0. 




Figure 6: Influence of roll angle on the aerodynamic coefficients, a = 5°, /5 = 7 = 
0. 
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Figure 7: Influence of fliglit patli angle on the aerodynamic coefficients, a = 5°, (3 
= = 0. 
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Figure 8: Vehicle rotational derivatives in ground effect in function of the angle of 
attack, p — (p — 0. 




Figure 9: Comparison of the dimensionless aerodynamic derivatives for h/c 
Data from Ref. 5 in horizontal axis, present results in vertical axis. 
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Figure 10: Trim calculation in horizontal flight at different heights. P 




Figure 11: Influence of tlie attitude on tlie trim in liorizontal fliglit. 




Figure 12: Trim analysis of a "truly banked" turn with ^ = 5° s 
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Figure 13: Root Locus. H = a) and b), H 7^ c) and d). 
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Figure 14: Time history of the perturbed longitudinal motion beginning from dif- 
ferent initial heights, (H 0). 
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Figure 15: Time history of the perturbed lateral motion beginning from different 
initial heights, (H 7^ 0). 




I I I I I I I I I I I I I I I I I I I I I I I I I 

2 4 6 8 10 



t (sec) 



44 



